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, x METHOD FOR REDUCING RINGING IN NUCLEAR MAGNETIC RESONANCE 

MEASUREMENTS 

§) BACKGROUND OF THE INVENTION 

Jv 

[0001] The present invention relates generally to a method for measuring nuclear 
magnetic resonance properties of an earth formation traversed by a borehole, and, more 
particularly, to a method for reducing ringing artifacts in a nuclear magnetic resonance 
measurement. 

[0002] Nuclear magnetic resonance (NMR) measurements taken in a borehole can 
provide different types of information about a geological formation. In the past, such 
measurements often were made after the borehole had been drilled. Today, however, it is 
possible to log NMR measurements while drilling (i.e., logging while drilling or LWD), thus 
saving time and providing valuable real-time information about the earth formation as drilling 
progresses. For example, such information can indicate the fractional volume of pore space, the 
fractional volume of mobile fluid, the total porosity of the formation, etc. 

[0003] NMR logging tools generally include one or more large permanent magnets 
or electromagnets for generating a static magnetic field, Bo, an antenna placed proximate the 
formation to be analyzed, and circuitry adapted to conduct a sequence of RF power pulses 
through the antenna to induce an RF magnetic field, Bi, in the formation. The circuitry also 
includes a receiver adapted to detect signals induced in the antenna as a result of the RF pulse 
sequence. The induced signals then can be measured and processed to provide the desired 
information about the properties of the formation. 

[0004] Typically, NMR logging tools are tuned to detect hydrogen resonance 
signals (e.g., from either water or hydrocarbons) because hydrogen nuclei are the most abundant 
and easily detectable. In general, measurements of NMR related phenomena of hydrogen nuclei 
in the earth formation are performed by allowing some time for the static magnetic field, B 0 > to 
polarize the hydrogen nuclei of water and hydrocarbons in a direction substantially parallel to B 0 , 
thus creating a nuclear magnetization. The direction of the nuclear magnetization can then be 
changed by applying a sequence of RF pulses to induce the RF field Bi. Commonly, the pulse 
sequence employed includes a first RF pulse (i.e., the excitation pulse) having a magnitude and 
duration selected to re-orient the nuclear magnetization by about 90 degrees from the orientation 
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attained as a result of B 0 (i.e., the initial transverse magnetization). After a selected time, a train 
of successive RF pulses is applied (i.e., inversion or refocusing pulses), each of which has a 
magnitude and a direction selected to re-orient the nuclear spin axes by about 1 80 degrees from 
their immediately previous orientations. The frequency of the RF field needed to re-orient the 
nuclear magnetization (i.e., the Larmor frequency) is related to the amplitude of the static 
magnetic field Bo by the gyromagnetic ratio y, which is unique to each isotope. 

[0005] After application of the initial RF pulse (i.e., after the nuclear magnetization 
is in the plane perpendicular to Bo), the nuclear magnetization begins to precess around Bo, 
producing a weak RF signal at the Larmor frequency which is detectable by the antenna Due to 
inhomogeneities in the magnetic field B 0 , the coherence between the individual spins eventually 
is lost and the nuclear magnetization decays rapidly. The inversion pulses re-create the lost 
magnetization (i.e., the coherence re-appears), producing signals that can be detected by the 
antenna. These signals, referred to as "spin echoes," generally are measured during the time 
interval between successive RF inversion pulses. The rate at which the spin echoes decay (i.e., 
the rate at which the nuclei irrevocably lose their alignment within the transverse plane) is 
referred to as the transverse relaxation rate. The time constant of this decay, referred to as the 
traverse relaxation time T2, is related to the chemical and physical properties of the earth 
formation. For example, hydrogen nuclei in viscous oils have relatively short relaxation times, 
whereas hydrogen nuclei in light oils have relatively long relaxation times. Similarly, hydrogen 
nuclei in free water typically have longer relaxation times than those in bound water (e.g., clay- 
bound water). 

[0006] To acquire the NMR data, several known pulse sequences are commonly 
employed. Such sequences include Carr-Purcell-like sequences, such as the Carr-Purcell- 
Meiboom-Gill ("CPMG") pulse sequence, which often is used for measuring transverse 
relaxation times. Advantages of the CPMG pulse sequence include compensation for pulse 
imperfections and inhomogeneities in both the B 0 and Bi magnetic fields, as is well known in the 
art. 

[0007] The Carr-Purcell-like sequences, however, do not compensate for all 
undesirable artifacts appearing in NMR measurements. For example, exciting the antenna with 
RF power pulses in the presence of a strong static magnetic field causes mechanical excitation of 
the antenna, which leads to generation of a spurious signal in the antenna called "ringing." The 
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ringing signal is unrelated to the NMR phenomena and typically has a large amplitude and slow 
rate of decay relative to the induced spin echo signals. Because the spin echoes typically are 
measured shortly after application of an RF pulse, the ringing signals often overwhelm the spin 
echo signals, resulting in loss of significant information regarding the formation properties. 

[0008] Several different techniques are known for reducing ringing. One technique 
has been to design the hardware to minimize the interaction between the electromagnetic fields 
and the materials of the hardware. An example of such hardware is described in U.S. Patent No. 
5,712,566 issued to Taicher et al. 

[0009] Another technique, known as "phase alternating pairs" or PAPs, includes 
creating a phase difference between the ringing signal and the spin echo signals, and summing or 
"stacking" multiple echo sequences to reduce the amplitude of the ringing signal in the final 
output For example, to compensate the ringing of the 180 degree pulses in a particular pulse 
sequence, the sequence may be repeated twice with opposing directions of the initial 90 degree 
pulse. By stacking the two sets of echo measurements, the ringing contribution of the 180 degree 
pulses can be substantially canceled. The ringing of the 90 degree pulse, on the other hand, is 
not canceled. However, because the 90 degree pulse is only applied once at the beginning of 
each sequence, its influence dies down relatively quickly and thus affects only early spin echo 
signals. Examples of ringing cancellation methods utilizing PAPs are described in U.S. Patent 
Nos. 5,596,274 issued to Sezginer and 5,023,551 issued to Kleinberg et al., and International 
Publication No. WO 98/43064 by Numar Corp. 

[0010] A necessary assumption underlying the phase alternated pair ringing 
cancellation technique is that the ringing does not change between the two measurement 
sequences. However, in borehole logging applications, the logging tool is moving continuously 
as the NMR measurements are being made. It has been found in such applications that 
movement and, in particular, bending of the tool changes the characteristic of the ringing signal 
between consecutive measurements. Because the two measurements forming a phase alternating 
pair are separated by a wait time of typically a few seconds, the orientation of the tool in the 
borehole during a first sequence bears no relationship to the orientation of the tool during the 
successive sequence. Thus, the phase alternating pair cancellation technique may not result in 
complete cancellation of the ringing signal in such applications. 
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[0011] Accordingly, it would be desirable to provide a new pulse sequence and a 
method of processing the spin echo signals resulting from the new pulse sequence to cancel the 
ringing. 

BRIEF SUMMARY OF THE INVENTION 

[0012] A method of reducing a ringing artifact in an NMR measurement is 
provided. In accordance with one aspect of the invention, the method comprises inducing a 
static magnetic field in a volume, inducing an RF magnetic field in the volume in accordance 
with a pulse sequence, and acquiring a sequence of signals generated in the volume in response 
to the pulse sequence. Each signal in the acquired sequence includes a spurious signal 
component and a spin echo component. At least a first signal and a second signal of the acquired 
sequence are combined, thereby generating a corrected signal having a reduced spurious signal 
component. 

[0013] In accordance with another aspect of the invention, a method of reducing a 
ringing signal generated while measuring a nuclear magnetic resonance property of an earth 
formation adjacent a borehole is provided. The method comprises inserting a logging tool into 
the borehole, applying a static magnetic field to polarize spins within a volume of the formation, 
and applying an RF magnetic field to the volume in accordance with a pulse sequence having a 
plurality of pulses. After each of the plurality of pulses, spin echoes induced in the volume are 
acquired, thereby forming a measurement set comprising a plurality of spin echo signals, each of 
which includes a ringing component. The spin echo signals within the measurement set are 
combined to reduce the ringing components, thereby generating a corrected measurement set. 

[0014] In accordance with yet another aspect of the invention, a method for 
determining an earth formation property from NMR measurements is provided. The method 
comprises acquiring a first measurement set of spin echo signals induced in the earth formation 
by a first RF pulse sequence, the spin echo signals including a spin echo component and a noise 
component. The spin echo signals are linearly combined to reduce the noise component, 
resulting in a first corrected measurement set. An NMR parameter is determined based on the 
first corrected measurement set, and the earth formation property then is derived from the NMR 
parameter. 
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[0015] The foregoing has outlined rather broadly the features and technical 
advantages of the present invention in order that the detailed description of the invention that 
follows may be better understood. Additional features and advantages of the invention will be 
described hereinafter which form the subject of the claims of the invention. It should be 
appreciated by those skilled in the art that the conception and specific embodiment disclosed 
may be readily utilized as a basis for modifying or designing other structures for carrying out the 
same purposes of the present invention. It should also be realized by those skilled in the art that 
such equivalent constructions do not depart from the spirit and scope of the invention as set forth 
in the appended claims. The novel features which are believed to be characteristic of the 
invention, both as to its organization and method of operation, together with further objects and 
advantages will be better understood from the following description when considered in 
connection with the accompanying figures. It is to be expressly understood, however, that each 
of the figures is provided for the purpose of illustration and description only and is not intended 
as a definition of the limits of the present invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0016] For a more complete understanding of the present invention, reference is 
now made to the following descriptions taken in conjunction with the accompanying drawing, in 
which: 

[0017] FIG. 1 is a diagram of a typical setup for well logging; 

[0018] FIG. 2 is a block diagram of an exemplary embodiment of circuitry for 
producing RF pulses and receiving the spin echoes; 

[0019] FIG. 3 generally illustrates a typical pulse sequence which includes an 
excitation pulse and multiple inversion pulses; 

[0020] FIG. 4 is a graph illustrating simulated data resulting from application of 
each of the CP, CPAP, CPMG, and anti-ringing pulse sequences; 

[0021] FIG. 5 is a graph illustrating actual data acquired as a result of application 
of each of the CP, CPAP, CPMG, and anti-ringing pulse sequences; 
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[0022] FIG. 6 is a graph illustrating sensitivity curves comparing porosity 
estimations derived from CPMG measurements and from two different applications of the anti- 
ringing measurement technique disclosed herein; and 

[0023] FIG. 7 is a graph comparing sensitivity curves for a porosity estimation 
derived from a traditional inversion technique using CPMG echoes, a echo reconstruction 
scheme, and a linear combination processing scheme. 

DETAILED DESCRIPTION OF THE INVENTION 

[0024] FIG. 1 illustrates a schematic of an exemplary nuclear magnetic resonance 
(NMR) well logging tool 10 for investigating one or more earth formations 12 traversed by a 
borehole 14. The logging tool 10 typically is suspended in the borehole 14 on an armored cable 
16, the length of which substantially determines the relative depth of the logging tool 10. The 
cable length is controlled by any suitable means, such as a drum and winch arrangement 1 8 at the 
surface. 

[0025] The NMR logging tool 10 can be any suitable downhole NMR logging 
device adapted for either wireline logging applications or logging-while-drilling (LWD) 
applications. As previously discussed, the logging tool 10 typically includes a permanent 
magnet, an electromagnet, or a magnet array for generating a static magnetic field B 0 in a volume 
of interest 19 in the formations, one or more RF antennae (e.g., solenoid antennae, loop 
antennae, saddle antennae, etc.), and circuitry configured to produce pulses of RF power to 
induce the RF magnetic field Bi in the formations and receive the spin echoes detected from the 
formations. 

[0026] A surface recording system 20 provides electrical power to the logging tool 
10, and signals detected by the tool 10 are returned to system 20 for recording and interpretation 
via the armored cable 16. Typically, the surface recording system 20 maintains a log of the 
detected spin echoes with respect to the depth of the logging tool 10. In the embodiment 
illustrated, output signals representative of depth are provided by a cable length measuring 
encoder 22. Alternatively, the surface recording system 20 may maintain a log of the detected 
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spin echoes with respect to time. Later, the time-based measurements may be correlated with a 
log of depth measurements such that depth-based measurements may be derived. 

[0027] FIG. 2 illustrates, in block diagram form, an exemplary embodiment of 
downhole circuitry associated with the logging tool 10 that is configured to produce the RF 
pulses and detect the spin echoes. In general, the circuitry generates an RF signal which is 
transmitted through an antenna to induce an RF magnetic field in the earth formations. Spin 
echo signals generated as a result of the RF magnetic field are detected by the antenna and either 
are stored or transmitted back to the earth's surface for logging by the surface recorder system 
20. One skilled in the art would appreciate that various embodiments of the circuitry may be 
used without departing from the scope of the invention. 

[0028] As shown in the exemplary embodiment of FIG. 2, the downhole circuitry 
includes a processor subsystem 210 having associated memory, timing circuitry, interfaces, and 
selected peripheral devices (not separately shown). The processor subsystem 210 is coupled 
with a telemetry circuitry 212, which communicates with the surface recording system 20. 

[0029] The pulse forming circuitry includes a variable frequency oscillator 214 
which, under control of the processor subsystem 210, generates an alternating RF signal at the 
desired frequency. The output of the oscillator 214 is coupled to a phase shifter 216 (which 
provides for control of pulse phases) and a modulator 218, both of which are under control of the 
processor subsystem 210 to produce the desired pulse phases of the RF field. The output of the 
modulator 218 is coupled, via a power amplifier 220, to an RF antenna 222. A Q-switch 224 
optionally can be provided to damp the RF antenna system to reduce antenna ringing. 

[0030] In the embodiment illustrated, the antenna 222 both transmits the RF pulses 
to induce the RF field in the formations and detects the echo signals resulting from application of 
the RF field. Thus, as shown in FIG. 2, the antenna 222 also is coupled with a receiver section 
through a duplexer or switch 226, the output of which is coupled with a receiver amplifier 228. 
During transmitting and damping modes, the switch 226 protects the receiver amplifier 228 
from the high power pulses which pass to the RF antenna 222. During the receiving mode, the 
duplexer 226 acts as a low impedance connection from the antenna 222 to the receiver amplifier 
228. The output of the receiver amplifier 228 is coupled with a dual phase-sensitive detector 
230, which also receives, as a reference, a signal derived from the oscillator signal. The output 
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of the detector 230 is coupled to an analog-to-digital converter 232, the output of which is a 
digital signal representative of the detected NMR signals. 

[0031] It should be understood that although the logging tool 10 is shown as an 
integral or unitary device in FIG. 1, it may alternatively comprise separate components and may 
be combinable with other logging tools. Further, while a wireline tool is illustrated in FIG. 1, 
alternative forms of physical support and a communicating link with the surface can be used, for 
example, in an LWD application. Still further, the digital signals representative of the detected 
spin echoes may be transmitted to the surface recording system 20 while the tool 10 is downhole. 
The transmitted signals may be either the received raw measurements or the relevant NMR 
parameters computed using downhole processing. Alternatively, the digital signals may be 
stored in memory by the processor subsystem 210 and later retrieved when the tool 10 has 
returned to the surface. 

[0032] Using the logging tool 10, several NMR parameters can be measured from 
which formation properties may be derived. For example, most NMR logging instruments are 
configured to measure the spin-lattice (longitudinal) relaxation times (7/) and/or spin-spin 
(transverse) relaxation times (T 2 ) of hydrogen nuclei. Such measurements may be acquired by 
first polarizing the nuclei by exposure to a static magnetic field B 0 and then applying an RF pulse 
(the excitation pulse) tuned at the Larmor frequency of the nuclei of interest and calibrated in 
length to achieve a 90 degree rotation of the spin magnetization. 

[0033] The following discussion is made with reference to a rotating frame in 
which the Bo field and the Bi field are constant. The rotating frame is defined such that its z-axis 
is parallel to Bo and its x-axis is parallel to Bi during the initial 90 degree pulse. Further, the 
following discussion is made with reference to FIG. 3, which generally illustrates a sequence of 
pulses that may be used to measure NMR properties. 

[0034] As shown in FIG. 3, there typically is a wait time T w (e.g., 2-3 seconds) 
before application of a 90 degree excitation pulse 234. Immediately after the excitation pulse 
234, the magnetization of the nuclei is aligned with they-axis of the rotating frame. The 
magnetization then decays rapidly, primarily due to unavoidable inhomogeneities of the 
magnetic field. This rapidly decaying signal is called the free induction decay, or FID, and is 
represented in FIG 3 by the signal 236. The FID typically is not detectable by well-logging 
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tools because of its fast decay and the ringing generated in the antenna as a result of the 90 
degree pulse. 

[0035] The signal can be re-created, however, by applying a 1 80 degree (i.e., 
inversion) pulse, such as the pulse 238 in FIG. 3, which causes the de-phased spins to regain 
phase coherence, thus forming a spin echo (e.g., signal 240). This process of echo formation by 
application of a 180 degree pulse can be repeated hundreds of times (i.e., a pulse sequence) to 
induce a corresponding sequence of spin echoes. In FIG. 3, the pulse sequence includes 
inversion pulses 238, 242, 244, 246, and so forth, and the corresponding spin echo sequence 
includes signals 240, 248, 250, 252, and so forth. Measurements of the spin echo signals 240, 
248, 250, and 252 are performed at intervals of t e , which typically is in the range of 0.5 to 1.5 
milliseconds. 

[0036] One well-known pulse sequence is the Carr-Purcell (CP) sequence, which 
can be symbolized by the following equation: 

+X!(-^)+X 2 (+3; 2 )+X3(-^)+ X4(+y<) ... (1) 

wherein the terms in parentheses denote the direction of echo formation, and the bold uppercase 
terms represent the direction in the rotating frame of the 180 degree inversion pulses. Thus, it 
can be seen from the foregoing equation that the transverse magnetization of the nuclei is 
inverted by each 180 degree pulse, resulting in spin echo signals formed in an alternating manner 
along the +y and -y direction of the rotating frame. Although the initial 90 degree pulse is not 
set forth in the foregoing equation, it is assumed that the pulse is parallel to the x-axis of the 
rotating frame and produces a free induction decay signal aligned with they-axis. 

[0037] The amplitudes of the spin echoes decay with a characteristic decay time T 2 , 
which theoretically could be measured with the CP pulse sequence. However, spin dynamics 
errors due to off-resonance effects and pulse imperfections resulting from unavoidable 
inhomogeneities in the static magnetic field B 0 and the RF magnetic field Bi accumulate during 
the CP pulse sequence, resulting in a rapid signal decay unrelated to T 2 . Consequently, it is not 
practically possible to measure T 2 using the CP pulse sequence. 

[0038] At least two alternative pulse sequences are known which can compensate 
for these otherwise accumulating effects. The first such sequence is a CP pulse sequence with 
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alternating pulses (designated for the purposes of this disclosure as CPAP), which is represented 
by the following equation (again, the excitation pulse is not represented): 

+X l (-y J yx 2 (+y 2 ^X 3 (-y 3 yX4^y4) . . . (2) 

In the CPAP sequence, the phase of the inversion pulses is reversed from one pulse to the next, 
thus providing compensation for spin dynamics errors. Accordingly, the otherwise rapid decay 
is substantially eliminated, enabling measurement of 7^. 

[0039] The second such alternative sequence is the well-known and commonly 
used Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence, which is represented by the following 
equation: 

^Y^y^^^i^Y^ ... (3) 

In the CPMG sequence represented by equation (3), a 90 degree phase shift (not shown) is 
inserted between the 90 degree excitation pulse (not shown) and the train of inversion pulses. 
Consequently, the inversion pulses all are either parallel or antiparallel to the transverse 
magnetization, and all spin echoes form along the same axis, thus compensating for otherwise 
accumulating spin dynamics errors. 

[0040] In the CP sequence represented by equation (1) above, consecutive spin 
echoes have alternating phases (i.e., alignment in the rotating frame), while all of the inversion 
pulses have the same phase. Thus, the relative phase between each spin echo and its preceding 
inversion pulse alternates between consecutive pulses. Consequently, because the phase of the 
ringing artifacts corresponds to the phase of the associated inversion pulse, the ringing artifacts 
can be canceled simply by combining spin echo signals within the sequence (assuming the 
ringing signal does not change between pulses which are separated by a time delay of the order 
of milliseconds). 

[0041] In the CPAP sequence represented by equation (2), both consecutive echoes 

and consecutive inversion pulses have alternating phases. Consequently, in the CPAP sequence, 

the relative phase between each spin echo and its preceding inversion pulse is constant. 

Similarly, in the CPMG sequence represented by equation (3), all inversion pulses have the same 

phase and produce spin echoes having the same phase, resulting in a constant relative phase 

between each spin echo and its preceding inversion pulse. Accordingly, for both CPAP and 
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CPMG sequences, ringing artifacts can be canceled only by applying two pulse sequences having 
excitation pulses with opposing phases and then combining or stacking the two resultant echo 
sequences (assuming that the ringing does not change between consecutive sequences which are 
separated by a wait time on the order of seconds). 

[0042] Although inter-sequence echoes resulting from the CP sequence can be 
linearly combined to cancel the ringing artifacts, an accurate measurement of T2 cannot be made 
because of the accumulating errors due to magnetic field inhomogeneities. And, although two 
pulse sequences resulting from CP AP and CPMG measurements can be stacked to cancel 
ringing, the assumption that ringing does not change between consecutive sequences may not 
hold true in certain logging applications (e.g., LWD applications) due to the relatively long wait 
time (e.g., seconds) between measurements. 

[0043] In particular, in many downhole NMR applications, the orientation of the 
logging tool is continuously changing as the measurements are being made. For example, in an 
LWD application, the tool is rotating while measurements are made and thus may be facing 
different formations having different characteristics resulting in different influences on the 
ringing signal. Further, the logging tool may bend during the logging operation, which also has 
been found to change the characteristics of the ringing signals. Accordingly, in practice, NMR 
measurements derived from stacked PAPs made using the CPAP and CPMG sequences still may 
include significant ringing artifacts. 

[0044] The various weaknesses inherent in utilization of the CP, CPAP, and CPMG 
sequences to cancel ringing in downhole NMR applications may be overcome by use of a new 
"anti-ringing" pulse sequence and inversion algorithm that minimize the possibility that tool 
orientation has changed significantly between measurements, while also providing compensation 
for accumulating errors due to magnetic field inhomogeneities. Thus, such an anti-ringing 
sequence preferably would offer the advantages of the CP sequence (i.e., the ability to linearly 
combine echoes closely spaced in time (e.g., milliseconds) to cancel ringing), which could be 
accomplished by inclusion of either phase alternating inversion pulses or phase alternating spin 
echoes. It would further be desirable if the ringing cancellation properties of the anti-ringing 
sequence could be combined with the advantages of the CPAP/CPMG sequences (i.e., 
compensation for spin dynamics errors). 
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[0045] An anti-ringing pulse sequence that provides both substantial ringing 
cancellation and spin dynamics compensation is represented by the following equation (in which 
the excitation pulse is not shown): 

+Xi(-yi)-Xa^^^ ... (4) 

This sequence is a repeating pattern of three inversion pulses in the +x direction and three 
inversion pulses in the -x direction. To compensate for spin dynamics effects, the sequence 
incorporates a phase alternation (i.e., -X 2 ) after the first echo (i.e., -yfi. The phase alternation is 
followed by two more inversion pulses in the -x direction and two inversion pulses in the +x 
direction. Accordingly, this sequence is referred to as a six-cycle sequence (i.e., the phase 
pattern of the inversion pulses repeats every six pulses). 

[0046] FIGS. 4 and 5 are graphs of data points resulting from application of the CP, 
CPAP, CPMG, and the anti-ringing sequence in equation (4). FIG. 4 represents the first two 
hundred spin echo signals derived from a simulation. FIG. 5 illustrates the first two hundred spin 
echo signals acquired in a laboratory experiment with an inter-echo spacing of 1 .2 msec. Both 
graphs illustrate and compare the capability of each pulse sequence to compensate for spin 
dynamics effects. For example, it can be seen from FIGS. 4 and 5 that the signal resulting from 
the CP sequence (254) decays quickly due to accumulating off-resonance and field 
inhomogeneity effects, thus rendering measurement of T2 practically impossible. In contrast, the 
new anti-ringing sequence (256) performs almost as well as the CPMG (258) and CPAP (260) 
sequences for two out of every three echo signals. As can be seen in FIGS. 4 and 5, every third 
echo signal induced by the anti-ringing sequence is weaker by a factor of approximately 0.6. 

[0047] With respect to the anti-ringing pulse sequence represented by equation (4), 
cancellation of inversion pulse ringing can be performed by linearly combining and averaging 
nearby echoes (i.e., echoes that are separated by a relatively short time t e ) in the same sequence. 
In one exemplary embodiment, the time separation t e between combined echoes is in the range of 
2-3 milliseconds. However, the time separation may be shorter or longer than 2-3 milliseconds 
as the goal of using closely spaced echoes simply is to reduce the possibility that the position of 
the logging tool has changed significantly between measurements and thus has altered the 
ringing component of the acquired echo signals. 
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[0048] For the foregoing sequence, several different linear combinations of the 
echo signals result in ringing cancellation. For example, at least the following four different 
linear combinations yield results adequate to obtain measurements of the T2 distribution: 
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where Ic represents the linear combination, and e represents the individual echo signals. 

[0049] Thus, by using the linear combinations, the anti-ringing pulse sequence 
represented by equation (4) can result in reduction of ringing contributions that arise from 
changes in antenna orientation (e.g., bending, direction) without the need to take a second 
measurement to form a phase alternating pair. This ringing reduction technique is particularly 
advantageous in an LWD application, because ringing cancellation calculations are performed by 
combining and averaging echoes that are separated in time by only a few milliseconds. Thus, the 
opportunity for the orientation of the logging tool to change significantly is minimized. For 
example, in a typical logging application, the logging tool moves at approximately 120 r.p.m 
Consequently, the tool will move only approximately 2 to 3 degrees between cancellation echoes 
that have a time separation of 2-3 msec. In contrast, if phase alternating pairs are formed to 
cancel ringing, the cancellation echoes will be acquired several seconds apart, practically 
assuring that the orientation of the logging tool (and thus the ringing contribution) will have 
changed significantly between measurements. 

[0050] Although the formation of phase alternating pairs is not necessary to obtain 
ringing reduction, PAPs may be used in combination with the ringing reduction opportunity 
afforded by the new pulse sequence. Indeed, combining the two ringing reduction techniques 
may be particularly desirable, because the new echo cancellation technique will reduce ringing 
contributions influenced by movement of the logging tool, while formation of PAPs will reduce 
those ringing contributions that are independent of tool orientation and/or that change slowly 
from echo to echo. 
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[0051] Other sequences that provide for both ringing reduction via linear 
combination of closely spaced echoes within a single sequence and spin dynamics compensation 
also are contemplated. For example, 

+Y l (+y I yY 2 (+y 2 hY 3 (+y 3 ) ... (5) 

and 

+ Yi(ty/>f Y 2 (+y2)-Y3(+y3) ... (6) 

both are three-cycle sequences which compensate for spin dynamics errors and allow for linear 
combination of echoes (note: the excitation pulses are not shown). For example, in equation (5), 
echo 1 may be combined with echo 2 and echo 3. In equation (6), echo 3 may be combined with 
echo 1 and echo 2 to cancel ringing In both sequences, phase alternation of at least one of the 
inversion pulses provides compensation for otherwise accumulating spin dynamics errors. 

[0052] Another exemplary sequence is set forth in equation (7): 

+Xi(-y 7 )+Y 2 B^^ ... (7) 

Equation (7) is a six-cycle sequence (note: the excitation pulse is not shown) which generally is 
similar to a CP sequence with Y inversion pulses (i.e., pulses 2 and 5) interspersed to 
compensate for spin dynamics errors. Useful linear combinations to cancel ringing include 
combinations of echo 1 and echo 3, echo 1 and echo 6, echo 3 and echo 4, and echo 4 and echo 

6. 

[0053] Yet another exemplary sequence is set forth in equation (8): 

+YiW 2 (+^^^ ... (8) 

Equation (8) is a six-cycle sequence (excitation pulse is not shown) which appears similar to a 
CPMG sequence and thus provides compensation for spin dynamics errors. Ringing cancellation 
is provided by the inclusion of inversion pulses Y 2 and Y3 which are directed along the -y-axis in 
the rotating frame. Thus, each of echoes y 2 and y$ may be linearly combined with each of echoes 
yi,y4,y5, andy 6 to cancel ringing. 
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[0054] As can be seen from equations (4)-(8) above, not all combinations of inter- 



sequence echoes will have an effect on ringing cancellation. That is, either the echoes should be 
oppositely phased or the associated inversion pulses should be oppositely phased if the 
combination is to have an effect on ringing cancellation Accordingly, to optimally use the 
ringing cancellation feature offered by the pulse sequences represented by equations (4)-(8) 
above, only those linear combinations of echoes that contribute to ringing cancellation may be 
used to derive the NMR characteristic of interest (e.g., T 2 distribution). It should be understood, 
however, that other linear combinations other than those that cancel ringing also may used in the 
inversion algorithms. In such a case, the information of interest still may be derived, although 
the results may not be as sensitive. 



porosity estimation derived from different types of measurements. The curve 300 represents the 
sensitivity curve using spin echo measurements obtained from a CPMG sequence. The curves 
302 and 304 represent sensitivity curves using linear combinations of spin echoes as described 
above. In deriving curve 302, linear combinations with and without ringing cancellation were 
used. In contrast, in deriving curve 304, only those linear combinations that contribute to ringing 
cancellation were allowed. As can be seen from FIG. 6, deriving the solution from CPMG 
measurements provides the best results, and restricting the solution to using only ringing 
canceling linear combinations provides only slightly worse results than not imposing the 
restriction. 

[0056] The effect of using linear combinations of echoes on a T 2 inversion 
algorithm is illustrated by the following discussion. When deriving T 2 , the inversion algorithm 
attempts to solve the following system of linear equations: 



[0055] For example, FIG. 6 illustrates sensitivity curves 300, 302, and 304 for a 



(9) 



where e represents the measured echoes, and a represents the sought T 2 distribution. 



[0057] The matrix M in equation (9) is defined as: 




(10) 
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where t e is the inter echo spacing, t w is the wait time, and T 2J for ally = 1 , ... , N T2 represents the 
N T 2 bin values logarithmically sampling the interesting T 2 range. 

[0058] The second factor set forth in equation (10) describes the polarization 
correction, where the function f(T 2j) = Tjj describes the longitudinal relaxation time 7/ 
corresponding to the respective T 2 . 

[0059] Equation (9) can be solved by a constrained minimization of: 

F(a)=\Ma-ef+Zp\\ 2 (11) 

where a regularization term with the regularization parameter X is added. 

[0060] To force the T 2 inversion algorithm to use the linear combinations of echoes 
instead of the raw echo measurements, a matrix X may be introduced, which describes the linear 
combinations of echoes and the spin dynamics correction factors. As a result, equation (9) may 
be rewritten as follows: 

Xe=XMa (12) 
and the cost function represented by equation (11) is rewritten as follows: 

F(5) = |X(M5-e)|| 2 +A||a|| 2 (13) 

with the following unconstrained solution: 

a = (M T X T XM + Xl)~ ] M T X T Xe (14) 

[0061] Based on this solution for the T 2 distribution, characteristics of the earth 
formation may be derived. For example, the porosity O of the earth formation can be computed 
as the sum of all elements of a . This computation is represented by the following equation: 

O = 13 = T(M T X T XM + Aiy 1 M r X T Xe = I T e (15) 
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where 2 = (1, 1, 1, ... ) is a row vector of "1" and the linear estimator I . The standard deviation 
of the thus computed porosity estimator is the standard deviation of the product of the individual 
echoes and the/,2 norm of / . 

[0062] The response s of the porosity estimator with respect to 7i then can be 
computed as follows: 

s — / T M (16) 

[0063] The foregoing example described a processing method which is optimal in 
the least squares sense. However, for some applications, it may be desirable to compute 
corrected echoes from the measured echoes that mimic CPMG echoes. If the corrected echoes 
mimic CPMG echoes, then commonly used inversion algorithms may be used to derive the 
desired NMR parameter, which could facilitate data delivery as well as interpretation of the data 

[0064] Corrected, or CPMG-like, echoes may be computed from the linear 
combinations of echoes. The linear combinations are computed from the measured echoes as 
follows: 

f=Xe (17) 

where e are the measure echoes. 

[0065] The formal solution for reconstructed echoes e R is: 

e R = (X T X)' l X r l = (X T Xy l X J e (18) 

However, this solution is problematic because the formation of linear combinations destroys 
some information (e.g., the ringing) which cannot be reconstructed from the linear combinations 
alone without making further assumptions. This solution also is problematic because it does not 
apply correction factors for spin dynamics. To apply the spin dynamics corrections without 
disturbing the ringing canceling properties of the linear combinations, the spin dynamics 
corrections must be applied to the linear combination rather than to the raw echo measurements. 
To accomplish this result, the linear combination formation matrix X in equation (18) may be 
replaced with the matrix X S d that includes the appropriate corrections. Thus, 
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(19) 



[0066] To reconstruct the echoes from the linear combinations alone, additional 
information must be supplied. Considering that the linear combinations cancel ringing, then the 
same linear combinations may be computed except that the echoes are subtracted instead of 
added. Thus, for example, the following combinations may be performed to reconstruct the 
echoes from the pulse sequence set forth in equation (4) above: 

(a) Tc x = - €j) V/ = 1, 4, 7, j = 3, 6, 9... 

(b) rc 2 =i(^,)V/ = 2,5,8,...,y=3,6,9,... 

(c) ^,=^-6^)^=3,6,9, ...,/=4,7,10,... 

(d) 7c 4 =±{e-e.)Vi = 3, 6,9,...,; =5,8,11,... 

[0067] These calculations result in a signal that primarily is the ringing signal with 
a very small NMR signal (caused by the T 2 decay of the NMR signal between the combined 
echoes). A reasonable constraint would be that these linear combinations equal zero, or: 

f = Xe=0 (20) 
where X describes the modified linear combinations. The inclusion of this constraint leads to: 



T 



XSD 1 

*» <21) 



and 

^=(x T X + X T x)~Vx OT e (22) 

[0068] Both the echo reconstruction scheme and the linear combination processing 
scheme set forth in equations (9)-(22) above have a poorer sensitivity than a traditional inversion 
computed from CPMG echoes. However, both of the linear combination schemes do have 
comparable sensitivity responses, indicating that the computation of reconstructed echoes from 
the linear combinations does not severely affect the inversion algorithm. This can be seen in 
FIG. 7, which compares the sensitivity of a porosity estimation out of CPMG echoes (curve 702), 
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the echo reconstruction scheme (curve 704), and the linear combination processing scheme 
(curve 706). 

[0069] It should be apparent from the foregoing description that many variations of 
pulse sequences can be used in the described ringing cancellation technique. Further, many 
variations of linear combinations of echoes within a sequence which reduce the ringing 
component in the spin echo signal also may be used, as well as different processing methods to 
derive the desired NMR parameter. Still further, the technique may also be suitable for 
canceling other noise or spurious components superimposed on the spin echo signal. It should 
also be understood that the technique is suited for use in both wireline applications and logging 
while drilling applications, as well as other types of NMR applications in which reduction of 
spurious signals is a concern. 

[0070] Although the present invention and its advantages have been described in 
detail, it should be understood that various changes, substitutions and alterations can be made 
herein without departing from the spirit and scope of the invention as defined by the appended 
claims. Moreover, the scope of the present application is not intended to be limited to the 
particular embodiments of the process, machine, manufacture, composition of matter, means, 
methods and steps described in the specification. As one of ordinary skill in the art will readily 
appreciate from the disclosure of the present invention, processes, machines, manufacture, 
compositions of matter, means, methods, or steps, presently existing or later to be developed that 
perform substantially the same function or achieve substantially the same result as the 
corresponding embodiments described herein may be utilized according to the present invention. 
Accordingly, the appended claims are intended to include within their scope such processes, 
machines, manufacture, compositions of matter, means, methods, or steps. 
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